Origin of the neutron skin thickness of 208 Pb in nuclear mean-field models 
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We study whether the neutron skin thickness Ar np of 208 Pb originates from the bulk or from the 
surface of the nucleon density distributions, according to the mean-field models of nuclear structure, 
and find that it depends on the stiffness of the nuclear symmetry energy. The bulk contribution to 
Ar np arises from an extended sharp radius of neutrons, whereas the surface contribution arises from 
different widths of the neutron and proton surfaces. Nuclear models where the symmetry energy 
is stiff, as typical relativistic models, predict a bulk contribution in Ar np of 208 Pb about twice as 
large as the surface contribution. In contrast, models with a soft symmetry energy like common 
nonrelativistic models predict that Ar np of 208 Pb is divided similarly into bulk and surface parts. 
Indeed, if the symmetry energy is supersoft, the surface contribution becomes dominant. We note 
that the linear correlation of Ar np of 208 Pb with the density derivative of the nuclear symmetry 
energy arises from the bulk part of Ar np . We also note that most models predict a mixed- type 
(between halo and skin) neutron distribution for 208 Pb. Although the halo-type limit is actually 
found in the models with a supersoft symmetry energy, the skin-type limit is not supported by 
any mean-field model. Finally, we compute parity-violating electron scattering in the conditions of 
the 208 Pb parity radius experiment (PREX) and obtain a pocket formula for the parity-violating 
asymmetry in terms of the parameters that characterize the shape of the 208 Pb nucleon densities. 
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I. INTRODUCTION 

The study of the size and shape of the density distri- 
butions of protons and neutrons in nuclei is a classic, yet 
always contemporary area of nuclear physics. The pro- 
ton densities of a host of nuclei are known quite well from 
the accurate nuclear charge densities measured in exper- 
iments involving the electromagnetic interaction like 
elastic electron scattering. In contrast, the neutron den- 
sities have been probed in fewer nuclei and are generally 
much less certain. 

The neutron distribution of 208 Pb, and its rms radius 
in particular, is nowadays attracting significant interest 
in both experiment and theory. Indeed, the neutron skin 
thickness, i.e., the neutron-proton rms radius difference 

Arv = (r 2 )^ - <r 2 >y 2 , (1) 

of this nucleus has close ties with the density-dependent 
nuclear symmetry energy and with the equation of state 
of neutron-rich matter. In nuclear models, Ar np of 208 Pb 
displays nearly linear correlations with the slope of the 
equation of state of neutron matter 0-3, with the den- 
sity derivative L of the symmetry energy [5l-lll| , and with 
the surface-symmetry energy of the finite nucleus 0] . At 
first sight, it may seem intriguing that a property of 
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the mean position of the surface of the nucleon densi- 
ties (Ar np ) is correlated with a purely bulk property of 
infinite nuclear matter (L). However, we have to keep in 
mind that Ar np depends on the surface-symmetry energy. 
This quantity reduces the bulk-symmetry energy due to 
the finite size of the nucleus. Assuming a local density 
approximation, we can correlate the surface-symmetry 
energy with the density slope L, which determines the de- 
parture of the symmetry energy from the bulk value. The 
correlation of Ar np with L then follows. Actually, these 
correlations can be derived almost analytically starting 
from the droplet model (DM) of Myers and Swiatecki 
[HI, [r| as we showed in Refs. [t| [l(|. By reason of 
its close connections with the nuclear symmetry energy, 
knowing accurately Ar np of 208 Pb can have important 
implications in diverse problems of nuclear structure and 
of heavy-ion reactions, in studies of atomic parity viola- 
tion, as well as in the description of neutron stars and in 
other areas of nuclear astrophysics (see, e.g., Refs. [l4T - 
l28l|). Since the charge radius of 208 Pb has been measured 
with extreme accuracy (r c h — 5.5013(7) fm [l|), the neu- 
tron rms radius of 208 Pb is the principal unknown piece 
of the puzzle. 

The lead parity radius experiment (PREX) [29] is a 
challenging experimental effort that aims to determine 
{t 2 )Ii ' 2 of 208 Pb almost model independently and to 1% 
accuracy by parity- violating electron scattering [2^, [3(| ■ 
This purely electroweak experiment has been run at the 
Jefferson Lab very recently, although results are not 
yet available. The parity-violating electron scattering 
is useful to measure neutron densities because in the 
low-momentum transfer regime the Z° boson couples 
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mainly to neutrons. For protons, this coupling is highly 
suppressed because of the value of the Weinberg angle. 
Therefore, from parity-violating electron scattering one 
can obtain the weak charge form factor and the closely 
related neutron form factor. From these data, the neu- 
tron rms radius can in principle be deduced [30]. This 
way of proceeding is similar to how the charge density is 
obtained from unpolarized electron scattering data [30|. 
The electroweak experiments get rid of the complexities 
of the hadronic interactions and the reaction mechanism 
does not have to be modeled. Thus, the analysis of the 
data can be both clean and model independent. There 
may be a certain model dependence, in the end, in hav- 
ing to use some neutron density shape to extract the 
neutron rms radius from the parity-violating asymmetry 
measured at a finite momentum transfer. 

To date, the existing constraints on neutron radii, 
skins, and neutron distributions of nuclei have mostly 
used strongly interacting hadronic probes. Unfortu- 
nately, the measurements of neutron distributions with 
hadronic probes are bound to have some model de- 
pendence because of the uncertainties associated with 
the strong force. Among the more frequent experimen- 
tal techniques we may quote nucleon elastic scattering 
[3ll [32( 1 , the inelastic excitation of the giant dipole and 
spin-di pole resonances [33l l34j , and experiments in exotic 
atoms |35l438| . Recent studies indicate that the pygmy 
dipole resonance may be another helpful tool to constrain 
neutron skins [HI, Ha]. 

The extraction of neutron radii and neutron skins from 
the experiment is intertwined with the dependence of 
these quantities on the shape of the neutron distribution 
[35l - l39( ] . The data typically do not indicate unambigu- 
osly, by themselves, if the difference between the periph- 
eral neutron and proton densities that gives rise to the 
neutron skin is caused by an extended bulk radius of the 
neutron density, by a modification of the width of the 
surface, or by some combination of both effects. In the 
present work we look for theoretical indications on this 
problem and study whether the origin of the neutron skin 
thickness of 208 Pb comes from the bulk or from the sur- 
face of the nucleon densities according to the mean-field 
models of nuclear structure. The answer turns out to 
be connected with the density dependence of the nuclear 
symmetry energy in the theory. 

We described in Ref. [4(| a procedure to discern bulk 
and surface contributions in the neutron skin thickness of 
nuclei. It can be applied to both theoretical and experi- 
mental nucleon densities as it only requires to know the 
equivalent sharp radius and surface width of these densi- 
ties, which one can obtain by fitting the actual densities 
with two-parameter Fermi (2pF) distributions. The 2pF 
shape is commonly used to characterize nuclear densities 
and nuclear potentials in both theoretical and experi- 
mental analyses. The doubly magic number of protons 
and neutrons in 208 Pb ensures that deformations do not 
influence the results and spherical density distributions 
describe the nuclear surface very well. We perform our 



calculations with several representative effective nuclear 
forces, namely, nonrelativistic interactions of the Skyrme 
and Gogny type and relativistic mean-field (RMF) in- 
teractions. The free parameters and coupling constants 
of these nuclear interactions have usually been adjusted 
to describe data that are well known empirically such 
as binding energies, charge radii, single-particle proper- 
ties, and several features of the nuclear equation of state. 
However, the same interactions predict widely different 
results for the size of neutron skin of 208 Pb and, as we will 
see, for its bulk or surface nature. We also study the halo 
or skin character [35-39] of the nucleon densities of 208 Pb 
in mean-field models. Finally, we perform calculations of 
parity-violating electron scattering on 208 Pb. We show 
that if 2pF nucleon densities are assumed, the parity- 
violating asymmetry as predicted by mean-field models 
can be approximated by a simple and analytical expres- 
sion in terms of the central radius and surface width of 
the neutron and proton density profiles. This suggests 
that an experiment such as PREX could allow to obtain 
some information about the neutron density profile of the 
208 Pb nucleus in addition to its neutron rms radius. 

The rest of the article proceeds as follows. In Sec.Qll we 
summarize the formalism to decompose the neutron skin 
thickness into bulk and surface components. The results 
obtained in the nuclear mean-field models are presented 
and discussed in Sec. lIIII A summary and the conclusions 
are given in Sec. IIVI 

II. FORMALISM 

The analysis of bulk and surface contributions to the 
neutron skin thickness of a nucleus requires proper defi- 
nitions of these quantities based on nuclear density dis- 
tributions. We presented in Ref. [5(| such a study, and 
we summarize only its basic points here. 

One can characterize the size of a nuclear density dis- 
tribution p(r) through several definitions of radii, and 
each definition may be more useful for a specific purpose 
(see Ref. [4l| for a thorough review). Among the most 
common radii, we have the central radius C: 

the equivalent sharp radius R: 

a r°° 

-7ri? 3 p(bulk) = 4vr / p(r)r 2 dr, (3) 
3 Jo 

i.e., the radius of a uniform sharp distribution whose den- 
sity equals the bulk value of the actual density and has 
the same number of particles; and the equivalent rms ra- 
dius Q: 

- 5 Q 2 = (r 2 ), (4) 

which describes a uniform sharp distribution with the 
same rms radius as the given density. These three radii 
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are related by expansion formulas |4l| : 



R{ 1 



Here, b is the surface width of the density profile: 



1 



P(0) 



ar 



(5) 



(6) 



which provides a measure of the extent of the surface 
of the density. Relations ([5]) usually converge quickly be- 
cause b/R is small in nuclei, especially in heavy-mass sys- 
tems. Nuclear density distributions have oscillations in 
the inner bulk region and a meaningful average is needed 
to determine the density values p(0) and ^(bulk) appear- 
ing in the above equations. This can be achieved by 
matching the original density with a 2pF distribution: 



p{r) = 



Po 



1 + cxp [(r - C)/a] 



(7) 



In 2pF functions the bulk density value corresponds very 
closely to the central density, and the latter coincides 
to high accuracy with the po parameter if exp (— C/a) is 
negligible. The surface width b and the diffuseness pa- 
rameter a of a 2pF function are related by b = (7r/\/3)a. 

As discussed in Ref . [4l[ , the equivalent sharp radius R 
is the quantity of basic geometric importance of the C, Q, 
and R radii. This is because a sharp distribution of radius 
R has the same volume integral as the density of the finite 
nucleus and differs from it only in the surface region. We 
illustrate this fact in Fig. [T] using as example the neutron 
density of 208 Pb of a mean-field calculation. We can see 
that the mean-field density is clearly overestimated in the 
whole nuclear interior by a sharp sphere of radius C. The 
equivalent rms radius Q fails also, by underestimating it. 
Only the equivalent sharp radius R is able to reproduce 
properly the bulk part of the original density profile of 
the nucleus. Therefore, R appears as the suitable radius 
to describe the size of the bulk of the nucleus. 

As the neutron skin thickness ([1]) is defined through 
rms radii, it can be expressed with Q: 



Ar np = \ - (Q n 



(8) 



Recalling from (0 that Q ~ R + |(6 2 /i?), we have a 



natural distinction in Ar np between bulk (oc R n 
and surface contributions. That is to say, 



with 
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independent of surface properties, and 
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' np 
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FIG. 1: (Color online) Comparison of sharp density profiles 
having radii C, R, and Q with the mean-field and 2pF density 
distributions for the neutron density of 20S Pb. The RMF 
interaction NL3 has been used in the mean-field calculation. 



of surface origin. The nucleus may develop a neutron 
skin by separation of the bulk radii R of neutrons and 
protons or by modification of the width b of the surfaces 
of the neutron and proton densities. In the general case, 
both effects are expected to contribute. We note that Eq. 
(fTTj) coincides with the expression of the surface width 
contribution to the neutron skin thickness provided by 
the DM of Myers and Swiatecki [H, EH if we set in Eq. 
CP R n = R P = roA 1 / 3 . 

The next-order correction to Eq. (fTTj) can be eas- 
ily evaluated for 2pF distributions (cf. Ref. [4l| for the 
higher-order corrections to the expansions (|5|)) and gives 



A „surr .corr 
np 



3 5 21/^^^\ 
5 2 20 Vig R 3j- 



(12) 



This quantity is usually very small — indeed, we neglected 
it in [40]. In the case of 208 Pb, we have found in all 
calculations with mean-field models that it is between 
—0.0025 fm and —0.004 fm, and thus can be obviated for 
most purposes. But because in the present work we deal 
with some detailed comparisons among the models, we 
have included (1121) in the numerical values shown for the 



r surf 
np 



It is to be mentioned that there is no universal method 
to do the parametrization of the neutron and proton den- 
sities with 2pF functions. A popular prescription is to use 
a x 2 minimization of the differences between the density 
to be reproduced and the 2pF profile, or of the differences 
between their logarithms. These methods may somewhat 
depend on conditions given during minimization (number 
of mesh points, limits, etc.). As in [ioj . we have preferred 
to extract the parameters of the 2pF profiles by imposing 
that they reproduce the same quadratic (r 2 ) and quartic 
(r 4 ) moments of the self-consistent mean-field densities, 
and the same number of nucleons. These conditions al- 
low us to determine in a unique way the equivalent 2pF 
densities and pay attention to a good reproduction of 
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the surface region of the original density because the lo- 
cal distributions of the quantities r 2 p(r) and r 4 p(r) are 
peaked at the peripheral region of the nucleus. An ex- 
ample of this type of fit is displayed in Fig. Q] by the 
dash-dotted line. It can be seen that the equivalent 2pF 
distribution nicely averages the quantal oscillations at 
the interior and reproduces accurately the behavior of 
the mean-field density at the surface. 



III. RESULTS 



Survey of model predictions and of data for 

Ar nv of 208 Pb 



We calculate our results with nonrelativistic models 
of Skyrmc type (SGII, Ska, SkM*, Sk-T4, Sk-T6, Sk- 
Rs, SkMP, SkSM*, MSk7, SLy4, HFB-8, HFB-17) and 
of Gogny type (D1S, DIN), as well as with several rela- 
tivistic models (NL1, NL-Z, NL-SH NL-RA1, NL3, TM1, 
NLC, G2, FSUGold, DD-ME2, NL3*). The original ref- 
erences to the different interactions can be found in the 
papers @,|42|for the Skyrme models, E3I for the Gogny 
models, and~|44T-|48| and Ref. [19] in HU for the RMF 
models. It may be mentioned that the recent force HFB- 
17 [42j achieves the lowest rms deviation with respect to 
experimental nuclear masses found to date in a mean- 
field approach. 

As is well known, nonrelativistic and relativistic mod- 
els differ in the stiffness of the symmetry energy. Note 
that by soft or stiff symmetry energy we mean that the 
symmetry energy increases slowly or rapidly as a func- 
tion of the nuclear density around the saturation point. 
Of course, the soft or stiff character can depend on the 
explored density region; for example, it is possible that a 
symmetry energy that is soft at nuclear densities becomes 
stiff at much higher densities (49|, or that a model with 
a stiff symmetry energy at normal density has a smaller 
symmetry energy at low densities The density de- 

pendence of the nuclear symmetry energy c sym (p) around 
saturation is frequently parametrized through the slope 
L of c sym (p) at the saturation density: 



L 



dc sym (p) 



dp 



(13) 



The pressure of pure neutron matter is directly propor- 
tional to L 51| and thus the L value has important im- 
plications for both neutron-rich nuclei and neutron stars. 
The symmetry energy of the Skyrme and Gogny forces 
analyzed in this work displays, as usual in the nonrela- 
tivistic models, from the very soft to the moderately stiff 
density dependence at nuclear densities (see Table Q] for 
the L parameter of the models) . On the contrary, the ma- 
jority of the relativistic parameter sets have a stiff or very 
stiff symmetry energy around saturation. The exception 
to the last statement in our case are the covariant pa- 
rameter sets FSUGold and DD-ME2 that have a milder 



TABLE I: Neutron skin thickness in 208 Pb calculated with the 
self-consistent densities of several nuclear mean-field models 
(Ar^p) and its partition into bulk and surface contributions 
defined in Sec. [TTJ as well as the relative weight of these bulk 
and surface parts. The models have been set out according 
to increasing Ar^p. The density slope L of the symmetry 
energy of the models is also listed. In order to help distinguish 
relativistic and nonrelativistic models, we have preceded the 
relativistic ones with an r in this table. 



Model 


(fm) 


(fm) 


(fm) 


bulk 

% 


surf 

% 


L 

(MeV) 


HFB-8 


n 

U 


1 1 K 
110 


0.031 


0.084 


27 


73 


14.8 


MSk7 


n 
u 


1 ID 


0.030 


0.086 


26 


74 


9.4 


D1S 


n 
u 


loO 


0.062 


0.073 


46 


54 


22.4 


SGII 


n 
u 


loO 


0.065 


0.071 


48 


52 


37.6 


DIN 


n 
u 


14Z 


0.070 


0.072 


49 


51 


31.9 


Sk-T6 


n 
u 


101 


0.067 


0.084 


44 


56 


30.9 


HFB-17 


n 
u 


1 K 1 
101 


0.067 


0.084 


44 


56 


36.3 


SLy4 


u 


1 1 : 1 
lOl 


0.086 


0.075 


53 


47 


46.0 


SkM* 





170 


0.093 


0.077 


55 


45 


45.8 


rDD-ME2 





193 


0.098 


0.095 


51 


49 


51.3 


SkSM* 





197 


0.116 


0.082 


58 


42 


65.5 


SkMP 





197 


0.123 


0.074 


62 


38 


70.3 


r FSUGold 





207 


0.105 


0.102 


51 


49 


60.5 


Ska 





211 


0.140 


0.071 


66 


34 


74.6 


Sk-Rs 





215 


0.146 


0.069 


68 


32 


85.7 


Sk-T4 





248 


0.163 


0.085 


66 


34 


94.1 


rG2 





257 


0.171 


0.086 


66 


34 


100.7 


rNLC 





263 


0.174 


0.089 


66 


34 


108.0 


r NL-SH 





266 


0.169 


0.097 


64 


36 


113.6 


rTMl 





271 


0.172 


0.098 


64 


36 


110.8 


r NL-RA1 





274 


0.179 


0.095 


65 


35 


115.4 


rNL3 





280 


0.185 


0.095 


66 


34 


118.5 


rNL3* 





288 


0.191 


0.097 


66 


34 


122.6 


rNL-Z 





307 


0.209 


0.098 


68 


32 


133.3 


rNLl 





321 


0.216 


0.105 


67 


33 


140.1 



symmetry energy than the typical RMF models. FSUG- 
old achieves this through an isoscalar-isovector nonlin- 
ear meson coupling [46[ and DD-ME2 because of having 
density-dependent meson-exchange couplings [47] ], 

In Table U we display the neutron skin thickness of 
208 Pb obtained from the self-consistent densities of the 
various interactions (denoted as Ar s np ). It is evident that 
the nuclear mean-field models predict a large window of 
values for this quantity. The nonrelativistic models with 
softer symmetry energies point toward a range of about 
0.1-0.17 fm. Most of the relativistic models, having a stiff 
symmetry energy, point toward larger neutron skins of 
0.25-0.3 fm. In between, the relativistic models DD-ME2 
and FSUGold predict a result close to 0.2 fm and the 
Skyrme interactions that have relatively stiffer symmetry 
energies fill in the range between 0.2 and 0.25 fm. 

Before proceeding, we would like to briefly survey some 
of the recent results deduced for Ar np in 208 Pb from ex- 
periment. For example, the recent analysis in Ref. [52] of 
the data measured in the antiprotonic 208 Pb atom |35ll36| 
gives Ar„ p = 0.16 ± (0.02) sta t ± (0.04) syst fm, including 
statistical and systematic errors. Another recent study 
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[53j of the antiprotonic data for the same nucleus leads 
to Ar np = 0.20 ± (0.04) oxp ± (0.05) th fm, where the theo- 
retical error is suggested from comparison of the models 
with the experimental charge density. These determi- 
nations are in consonance with the average value of the 
hadron scattering data for the neutron skin thickness of 
208 Pb, namely, Ar np ~ 0.165 ± 0.025 fm (taken from 
the compilation of hadron scattering data in Fig. 3 of 
Ref. [36j]). We may also mention that the constraints on 
the nuclear symmetry energy derived from isospin diffu- 
sion in heavy ion collisions of neutron-rich nuclei suggest 
Ar np = 0.22 ± 0.04 fm in 208 Pb H. Following Ref. 
(55l| . the same type of constraints exclude Ar np values 
in 208 Pb less than 0.15 fm. A recent prediction based 
on measurements of the pygmy dipole resonance in 68 Ni 
and 132 Sn gives Ar np = 0.194 ± 0.024 fm in 208 Pb [27[. 
Finally, we quote the new value Ar np = 0.211±g;g|| fm 
determined in [5r| from proton elastic scattering. Thus, 
in view of the empirical information for the central value 
of Ar np and in view of the Arf; p values predicted by the 
theoretical models in Table U it may be said that those 
interactions with a soft (but not very soft) symmetry en- 
ergy, for example, HFB-17, SLy4, SkM*, DD-ME2, or 
FSUGold, agree better with the determinations from ex- 
periment. Nevertheless, the uncertainties in the available 
information for Ar np are rather large and one cannot rule 
out the predictions by other interactions. If the PREX 
experiment [1^, [3(| achieves the purported goal of accu- 
rately measuring the neutron rms radius of 208 Pb, it will 
allow to pin down more strictly the constraints on the 
neutron skin thickness of the mean-field models. 



B. Bulk and surface contributions to Ar„ p of 208 Pb 
in nuclear models and the symmetry energy 

We next discuss the results for the division of the 
neutron skin thickness of 208 Pb into bulk (ArJ^ lk ) and 
surface (Ar! 1 " rf ) contributions in the nuclear mean-field 
models, following Sec. |IIJ We display this information in 
Table |U It may be noticed that the value of Ar^ lp lk plus 
Ar™ rf (quantities obtained from Eqs. (|TU)) -([l"2" j) ) agrees 
excellently with Ar^ p (neutron skin thickness obtained 
from the self-consistent densities). One finds that the 
bulk contribution ArJ^ lk to the neutron skin of 208 Pb 
varies in a window from about 0.03 fm to 0.22 fm. The 
surface contribution Ar™ rf is comprised approximately 
between 0.07 fm and 0.085 fm in the nonrelativistic 
forces, and between 0.085 fm and 0.105 fm in the rela- 
tivistic ones. Thus, whereas the bulk contribution to the 
neutron skin thickness of 208 Pb changes largely among 
the different mean-field models, the surface contribution 
remains confined to a narrower band of values. 

Table U shows that the size of the neutron skin thick- 
ness of 208 Pb is divided into bulk and surface contribu- 
tions in almost equal parts in the nuclear interactions 
that have soft symmetry energies (say, L ~ 20-60). This 
is the case of multiple nonrelativistic interactions and 



TABLE II: Equivalent sharp radius and surface width of the 
2pF neutron and proton density distributions of 208 Pb in 
mean-field models. Units are fm. 



Model 


Rn 


R P 


b n 


bp 


HFB-8 


6.822 


6.782 


0.991 


0.819 


MSk7 


6.847 


6.808 


0.980 


0.801 


D1S 


6.830 


6.751 


0.994 


0.846 


SGII 


6.890 


6.806 


0.971 


0.821 


DIN 


6.845 


6.755 


0.979 


0.828 


Sk-T6 


6.862 


6.775 


0.994 


0.820 


HFB-17 


6.883 


6.797 


0.996 


0.821 


SLy4 
SkM* 


6.902 


6.790 


1.007 


0.852 


6.907 


6.786 


1.007 


0.847 


DD-ME2 


6.926 


6.800 


1.026 


0.829 


SkSM* 


6.955 


6.805 


0.970 


0.790 


SkMP 


6.943 


6.784 


0.997 


0.839 


FSUGold 


6.971 


6.836 


1.024 


0.808 


Ska 


6.970 


6.789 


0.998 


0.844 


Sk-Rs 


6.950 


6.762 


0.962 


0.806 


Sk-T4 


6.991 


6.780 


1.008 


0.823 


G2 


7.037 


6.817 


1.012 


0.824 


NLC 


7.087 


6.863 


1.016 


0.820 


NL-SH 


7.039 


6.821 


0.989 


0.772 


TM1 


7.085 


6.862 


1.005 


0.787 


NL-RA1 


7.065 


6.834 


1.008 


0.797 


NL3 


7.060 


6.821 


1.017 


0.807 


NL3* 


7.052 


6.806 


1.026 


0.814 


NL-Z 


7.134 


6.865 


1.058 


0.847 


NL1 


7.100 


6.822 


1.065 


0.840 



of the covariant DD-ME2 and FSUGold parameter sets. 
When the symmetry energy becomes softer, the bulk part 
tends to be smaller. Indeed, we see that in the mod- 
els that have a very soft symmetry energy (L < 20), 
which we may call "supersoft" [57j, the surface contri- 
bution takes over and it is responsible for the largest 
part (~ 75%) of Ar np of 208 Pb. At variance with this 
situation, in the models with stiffer symmetry energies 
(L > 75) about two thirds of Ar np of 208 Pb come from 
the bulk contribution, as seen in the Skyrme forces of 
stiffer symmetry energy and in all of the relativistic forces 
that have a conventional isovector channel (G2, TM1, 
NL3, etc.). We therefore note that in a heavy neutron- 
rich nucleus with a sizable neutron skin such as 208 Pb, 
the nuclear interactions with a soft symmetry energy pre- 
dict that the contribution to Ar np produced by differing 
widths of the surfaces of the neutron and proton den- 
sities (b n b p ) is similar to, or even larger than, the 
effect from differing extensions of the bulk of the nucleon 
densities (R n 7^ R p ). On the contrary, the nuclear inter- 
actions with a stiff symmetry energy favor a dominant 
bulk nature of the neutron skin of 208 Pb, and then the 
largest part of Ar np is caused by R n ^ R p . We collect 
in Table |TT] the found equivalent sharp radii R n and R p 
and surface widths b n and b p of the densities of 208 Pb in 
the present mean-field models. 

As we have had the opportunity to see, the neutron 
skin thickness of a heavy nucleus is strongly influenced 
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by the density derivative L of the symmetry energy. In- 
deed, one easily suspects from Table|I]that Ar np of 208 Pb 
is almost linearly correlated with L in the nuclear mean- 
field models, which Fig. [2] confirms for the present in- 
teractions. The correlation of the neutron skin thickness 
of 208 Pb with L has been amply discussed in the litera- 
ture [5T4TT[. as it implies that an accurate measurement 
of the former observable could allow to tightly constrain 
the density dependence of the nuclear symmetry energy. 
In particular, we studied the aforementioned correlation 
in Ref. [H where it is shown that the expression of the 
neutron skin thickness in the DM of Myers and Swiatecki 
[l~2l HH can be recast to leading order in terms of the L 
parameter. To do that, we use the fact that in all mean- 
field models the symmetry energy coefficient computed 
at p 0.10 fm -3 is approximately equal to the DM sym- 
metry energy coefficient in 208 Pb which includes bulk- 
and surface-symmetry contributions [9(. In the standard 
DM, where the surface widths of the neutron and proton 
densities are taken to be the same [H, [l3[ , the neutron 
skin thickness is governed by the ratio between the bulk- 
symmetry energy at saturation J = c s ™U) and the sur- 
face stiffness coefficient Q of the DM [9|, [l0[ (the latter is 
not to be confused with the equivalent rms radius Q of 
Eq. ((4])). The DM coefficient Q measures the resistance 
of the nucleus against the separation of the neutron and 
proton surfaces to form a neutron skin. We have shown 
!). 10 in mean-field models that the DM formula for the 
neutron skin thickness in the case where one assumes 
b n = bp, undershoots the corresponding values computed 
by the semiclassical extended Thomas-Fermi method in 
finite nuclei and, therefore, a nonvanishing surface contri- 
bution is needed to describe more accurately the mean- 
field results. However, this surface contribution has a 
more involved dependence on the parameters of the in- 
teraction and does not show a definite correlation with 
the J/Q ratio (see Fig. 4 of Ref. OH). 

Now, we wondered to which degree the correlation with 
L of the neutron skin thickness of 208 Pb holds in its bulk 
and surface parts extracted from actual mean-field densi- 
ties. From our discussion of the indications provided by 
the DM, we can expect this correlation to be strong in 
the bulk part and weak in the surface part. Indeed, the 
plots of ArJ^ lk and Ar™ rf against L in Fig. [5] show that 
the bulk part displays the same high correlation with L 
as the total neutron skin (the linear correlation factor is 
of 0.99 in both cases), whereas the surface part exhibits 
a mostly flat trend with L. The linear fits in Fig. [2]of the 
neutron skin thickness of 208 Pb and of its bulk part have 
also quite similar slopes. One thus concludes that the lin- 
ear correlation of Ar np of 208 Pb with the density content 
of the nuclear symmetry energy arises mainly from the 
bulk part of Ar np . In other words, the correlation arises 
from the change induced by the density dependence of 
the symmetry energy in the equivalent sharp radii of the 
nucleon density distributions of 208 Pb rather than from 
the change of the width of the surface of the nucleon 
densities. 
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FIG. 2: (Color online) Correlation of the neutron skin thick- 
ness of 208 Pb and of its bulk and surface parts with the density 
derivative L of the nuclear symmetry energy. 



The value of about 0.1 fm that the surface contri- 
bution to Ar np takes in 208 Pb can be understood as 
follows starting from Eq. (|TT|) . Taking into account 
that in 2pF distributions fitted to mean-field densities 
R n ~ R p ~ 1.16A 1 / 3 fm and b n + b p ~ 1.8 fm (see Table 



IrE]) , Eq. (fTT]) can be approximated as 



Ar surt 



3A- 1 / 3 (6„-6 p ). 



(14) 



Given that b n — b p ~ 0.2 fm for 208 Pb on the average in 
mean-field models (see Table ITT|) . one finds Ar™ rf - 0.1 
fm, rather independently of the model used to compute 
it. It is interesting why the range of variation of b n with 
respect to b p is not larger in nuclear models, in view of 
the fact that R n — R v can take more different values. 
As discussed in Ref. [3fJ, this constraint is imposed on 
the models most likely by the mass fits. For example, a 
model having nucleon densities with very small or very 
large surface widths (i.e., very sharp or very extended 
surfaces) would provoke a large change in the surface 
energy of the nucleus, but that hardly would be successful 
to reproduce the known nuclear masses. 



C. Discussion of the shape of the neutron density 
profiles 

The use of 2pF functions in order to represent the nu- 
clear densities by approximate distributions is also quite 
common in the experimental investigations. The parame- 
ters of the proton 2pF distribution can be assumed known 
in experiment, by unfolding from the accurately mea- 
sured charge density [40|. However, the shape of the neu- 
tron density is more uncertain, and even if the neutron 
rms radius is determined, it can correspond to different 
shapes of the neutron density. Actually, the shape of the 
neutron density is a significant question in the extrac- 
tion of nuclear information from experiments in exotic 
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TABLE III: Central radius and surface diffuseness of the 2pF 
neutron and proton density distributions of 208 Pb in mean- 
field models. Units are fm. 



Model 


C n 


Cp 


a n 


a v 


Cn Cp 


a„ — a 


HFB-8 


6.679 


6.683 


0.546 


0.451 


-0.004 


0.095 


MSk7 


6.707 


6.714 


0.540 


0.442 


-0.007 


0.099 


D1S 


6.687 


6.649 


0.546 


0.464 


0.038 


0.082 


SGII 


6.753 


6.707 


0.536 


0.453 


0.046 


0.083 


DIN 


6.705 


6.657 


0.537 


0.453 


0.048 


0.084 


Sk-T6 


6.718 


6.676 


0.548 


0.452 


0.042 


0.096 


HFB-17 


6.739 


6.697 


0.549 


0.453 


0.042 


0.096 


SLy4 


6.755 


6.683 


0.555 


0.470 


0.072 


0.085 


SkM* 


6.760 


6.681 


0.555 


0.467 


0.079 


0.088 


DD-ME2 


6.774 


6.699 


0.566 


0.457 


0.075 


0.109 


SkSM* 


6.819 


6.713 


0.535 


0.436 


0.106 


0.099 


SkMP 


6.799 


6.680 


0.550 


0.463 


0.119 


0.087 


FSUGold 


6.821 


6.740 


0.564 


0.446 


0.081 


0.118 


Ska 


6.827 


6.684 


0.550 


0.465 


0.143 


0.085 


Sk-Rs 


6.817 


6.665 


0.530 


0.444 


0.152 


0.086 


Sk-T4 


6.846 


6.681 


0.555 


0.453 


0.165 


0.102 


G2 


6.891 


6.717 


0.558 


0.454 


0.174 


0.104 


NLC 


6.941 


6.765 


0.560 


0.452 


0.176 


0.108 


NL-SH 


6.900 


6.733 


0.546 


0.426 


0.167 


0.120 


TM1 


6.942 


6.772 


0.554 


0.434 


0.170 


0.120 


NL-RA1 


6.921 


6.741 


0.556 


0.440 


0.180 


0.116 


NL3 


6.914 


6.726 


0.560 


0.445 


0.188 


0.115 


NL3* 


6.903 


6.709 


0.566 


0.449 


0.194 


0.117 


NL-Z 


6.977 


6.761 


0.584 


0.467 


0.216 


0.117 


NL1 


6.940 


6.718 


0.587 


0.463 


0.222 


0.124 



ato ms [35l - l38j and from parity- violating electron scatter- 
ing [39(. To handle the possible differences in the shape 
of the neutron density when analyzing the experimental 
data, the so-called "halo" and "skin" forms are frequently 
used [35 39] . In the "halo- type" distribution the nucleon 
2pF shapes have C n = C p and a n > a p , whereas in the 
"skin- type" distribution they have a n = a p and C n > C p . 
To complete our study, we believe worth discussing the 
predictions of the theoretical models for the parameters 
of the 2pF shapes in 208 Pb. 

We compile in Table lllTl the central radii C n and C p and 
the diffuseness parameters a n and a p of the 2pF nucleon 
density profiles of 208 Pb obtained from the mean-field 
models of Table HI We see that C n of neutrons spans a 
range of approximately 6.7-6.85 fm in the nonrelativistic 
interactions and that it is of approximately 6.8-7 fm in 
the relativistic parameter sets. In the case of the proton 
density distribution, the value of C p is smallest (~6.65 
fm) in the two Gogny forces, it is about 6.67-6.71 fm in 
the Skyrme forces, and it is in a range of 6.7-6.77 fm 
in the RMF models. Then, we note that not only C n of 
neutrons but also C p of protons is generally smaller in the 
nonrelativistic forces than in the relativistic forces. The 
total spread in C p among the models (about 0.12 fm) is, 
though, less than half the spread found in C n (about 0.3 
fm). Indeed, the accurately known charge radius of 208 Pb 
is an observable that usually enters the fitting protocol 
of the effective nuclear interactions. 



If we inspect the results for the surface diffuseness of 
the density profiles of 208 Pb in Table Hill we see that a n of 
neutrons lies in a window of 0.53-0.59 fm (with the ma- 
jority of the models having a n between 0.545 and 0.565 
fm). The nonrelativistic interactions favor a„ < 0.555 
fm, whereas the RMF sets favor a n > 0.555 fm. This 
indicates that the fall-off of the neutron density of 208 Pb 
at the surface is in general faster in the interactions with 
a soft symmetry energy than in the interactions with 
a stiff symmetry energy. The surface diffuseness a p of 
the proton density spans in either the nonrelativistic or 
the relativistic models almost the same window of values 
(0.43-0.47 fm; with the majority of the models having a p 
between 0.445 and 0.465 fm). This fact is in contrast to 
the other 2pF parameters discussed so far. Actually, the 
a p value of the proton density can be definitely larger in 
some nonrelativistic forces than in some relativistic forces 
(for example, in the case of SkM* and NL3). One finds 
that the total spread of a n and a p within the analyzed 
models is quite similar: about 0.05 fm in both a n and a p . 
This spread corresponds roughly to a 10% variation com- 
pared to the mean values of a n and a p . It is remarkable 
that while among the models C n has a significantly larger 
spread than C p , the surface diffuseness a n of the neutron 
density has essentially the same small spread as the sur- 
face diffuseness a p of the proton density. As we have 
discussed at the end of Sec. IIII B[ this is likely imposed 
by the nuclear mass fits. It means that our ignorance 
about the neutron distribution in 208 Pb does not seem 
to produce in the mean-field models a larger uncertainty 
for a n of neutrons than for a p of protons, and that most 
of the uncertainty goes to the value of C n . 

The difference C n — C p of the central radii of the nu- 
cleon densities of 208 Pb turns out to range approximately 
between 0. and 0.2 fm. It is smaller for soft symmetry 
energies and larger for stiff symmetry energies. We real- 
ize that the limiting situation of a halo-type distribution 
where the nucleon densities of 208 Pb have C n = C p and 
a n > a p is actually attained in the nuclear mean-field 
models with a very soft symmetry energy (like in HFB-8 
or MSk7 where C n — C p even is slightly negative). The 
difference a n — a p of the neutron and proton surface dif- 
fuseness in 208 Pb is comprised between nearly 0.08 and 
0.1 fm in the nonrelativistic forces and between nearly 
0.1 and 0.12 fm in the RMF forces. This implies that no 
interaction predicts a n — a p of 208 Pb as close to vanishing 
as C n — C p is in some forces. Thus, the limiting situation 
where the nucleon densities in 208 Pb would have a n = a p 
and C n > C p is not found in the nuclear mean-field mod- 
els. Indeed, we observe in Table Mil that if C n — C p be- 
comes larger in the models, also a n — a p tends overall to 
become larger. In order to help visualize graphically the 
change in the mean-field nucleon densities of 208 Pb from 
having a nearly vanishing C n — C p or a large C n — C p , 
we have plotted in Fig. [3] the example of the densities 
of the MSk7 and NL3 interactions. On the one hand, 
we see that both models MSk7 and NL3 predict basi- 
cally the same proton density as expected. On the other 
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parity- violating asymmetry 



A L r = 





do_ 






da + 


da_ 


~dtt 


+ dfi 



(15) 



where da±/dQ is the elastic electron- nucleus cross sec- 
tion. The plus (minus) sign accounts for the fact that 
electrons with a positive (negative) helicity state scatter 
from different potentials (V±(r) = Vboulomb (r) ± Kreak (r) 
for ultra-relativistic electrons). Assuming for simplicity 
the plane wave Born approximation (PWBA) and ne- 
glecting nucleon form factors, the parity-violating asym- 
metry at momentum transfer q can be written as [301 ] 



FIG. 3: (Color online) Comparison of the nucleon densities 
predicted in 208 Pb by the mean-field models MSk7 (C n -C p ~ 
fm) and NL3 (C„ - C p w 0.2 fm). 



jPWBA 
A LR 



G F q 



2 r 



C p in MSk7 



C p in NL3 can be appreciated in the higher 



hand, the difference between having C, 
and C n > 

bulk and the faster fall-off at the surface of the neutron 
density of MSk7 compared with NL3. 

In summary, we conclude that in 208 Pb the nuclear 
mean-field models favor the halo-type distribution with 
C n ~ Cp and a n > a p if they have a very soft ( "super- 
soft") symmetry energy, they favor a mixed-type distri- 
bution if they have mild symmetry energies, and a situa- 
tion where C n is clearly larger than C p if the symmetry 
energy is stiff, but that the pure skin-type distribution 
where a n — a p = in 208 Pb is not supported (not even 
o>n — ftp ~ 0) by the mean-field models. Although the ex- 
perimental evidence available to date on the neutron skin 
thickness of 208 Pb is compatible with the ranges of the 
C n — Cp and a n — a p parameters considered in our study, 
it is not to be excluded that the description of a precision 
measurement in 208 Pb may need of nucleon densities with 



C n — C p or a n — a p values not fitting Table UTT1 However, a 
sizable deviation (such — 0) could mean that 

there is some missing physics in the isospin channel of 
present mean-field interactions, because once these inter- 
actions are calibrated to reproduce the observed binding 
energies and charge radii of nuclei they typically lead to 
the ranges of Table IHII 



D. Application to parity-violating electron 
scattering 

Parity-violating electron scattering is expected to be 
able to accurately determine the neutron density in a 
nucleus since the Z° boson couples mainly to neutrons 
[H Ho)]. Specifically, the PREX experiment HH aims 
to provide a clean measurement of the neutron radius 
of 208 Pb. In this type of experiments one measures the 



4 sin Q w 



F n {q)-F p (q) 



F P (q) 



(16) 



where sin 2 6w m 0.23 for the Weinberg angle and F n (q) 
and F p (q) are the form factors of the point neutron and 
proton densities. Because F p {q) is known from elastic 
electron scattering, it is clear from ([T6|) that the largest 
uncertainty to compute Alr comes from our lack of 
knowledge of the distribution of neutrons inside the nu- 
cleus. PREX intends to measure Alr in 208 Pb with a 
3% error (or smaller). This accuracy is thought to be 
enough to determine the neutron rms radius with a 1% 
error [H,^. 

To compute the parity-violating asymmetry we essen- 
tially follow the procedure described in Ref. [30]. For 
realistic results, we perform the exact phase shift anal- 
ysis of the Dirac equation for electrons moving in the 
potentials V±(r) [58|]. This method corresponds to the 
distorted wave Born approximation (DWBA). The main 
input needed for solving this problem are the charge and 
weak distributions. To calculate the charge distribution, 
we fold the mean-field proton and neutron point-like den- 
sities with the electromagnetic form factors provided in 
[59l ] . For the weak distribution, we fold the nucleon point- 
like densities with the electric form factors reported in 
[30I ] for the coupling of a Z° to the protons and neu- 
trons. We neglect the strange form factor contributions 
to the weak density [3(| ■ Because the experimental anal- 
ysis may involve parametrized densities, in our study we 
use the 2pF functions extracted from the self-consistent 
densities of the various models. The difference between 
Alr calculated in 208 Pb with the 2pF densities and with 
the self-consistent densities is anyway marginal at most. 
In Fig. SI AYr BA obtained with the Fermi distributions 
listed in Table IIIII is plotted against the values of C n — C p 
(lower panel) and a n — a p (upper panel). To simulate the 
kinematics of the PREX experiment [29[ , we set the elec- 
tron beam energy to 1 GeV and the scattering angle to 
5°, which corresponds to a momentum transfer in the 
laboratory frame of q — 0.44 fm -1 . 

First, one can see from Fig.@]that the mean-field calcu- 
lations constrain in a rather narrow window the value of 
the parity-violating asymmetry in 208 Pb. The increasing 
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FIG. 4: (Color online) Parity-violating asymmetry for 1 GeV 
electrons at 5° scattering angle calculated from the 2pF neu- 
tron and proton density distributions of 208 Pb in nuclear 
mean-field models. 



trend of Alr with decreasing C n — C p indicates that Alr 
is larger when the symmetry energy is softer. Note that 
a large value A LR : rs 7 x 10~ 7 (at 1 GeV and 5°) would be 
in support of a more surface than bulk origin of the neu- 
tron skin thickness of 208 Pb and of the halo- type density 
distribution for this nucleus. Second, AYr BA displays 
in good approximation a linear correlation with C n — C p 
(r = 0.978), while the correlation with a n — a p is not 
remarkable. Nevertheless, we have found a very good de- 
scription of AY^ BA of the mean-field models — well below 
the 3% limit of accuracy of the PREX experiment — by 
means of a fit in C n — C p and a n — a p (red crosses in 
Fig. g]): 



A fit 



vary C n and a n in (|17p to match the measured value. 
According to the predictions of the models in Table IIII1 
it would be reasonable to restrict this search to windows 
of about 0-0.22 fm for C n - C p and 0.08-0.125 fm for 
a« — o, p . Therefore, the result of a measurement of the 
parity-violating asymmetry together with Eq. (1171) (or 
Fig. 3]) would allow not only to estimate the neutron rms 
radius of 208 Pb but also to obtain some insight about 
the neutron density profile in this nucleus. This assumes 
that the experimental value for A^r will fall in, or at least 
will be not far from, the region allowed by the mean-field 
calculations at the same kinematics. 



IV. SUMMARY 

We have investigated using Skyrme, Gogny, and rel- 
ativistic mean-field models of nuclear structure whether 
the difference between the peripheral neutron and proton 
densities that gives rise to the neutron skin thickness of 
208 Pb is due to an enlarged bulk radius of neutrons with 
respect to that of protons or, rather, to the difference 
between the widths of the neutron and proton surfaces. 
The decomposition of the neutron skin thickness in bulk 
and surface components has been obtained through two- 
parameter Fermi distributions fitted to the self-consistent 
nucleon densities of the models. 

Nuclear models that correspond to a soft symmetry en- 
ergy, like various nonrelativistic mean-field models, favor 
the situation where the size of the neutron skin thickness 
of 208 Pb is divided similarly into bulk and surface com- 
ponents. If the symmetry energy of the model is "super- 
soft", the surface part even becomes dominant. Instead, 
nuclear models that correspond to a stiff symmetry en- 
ergy, like most of the relativistic models, predict a bulk 
component about twice as large as the surface compo- 
nent. We have found that the size of the surface compo- 
nent changes little among the various nuclear mean-field 
models and that the known linear correlation of Ar np 

208t 



[a + p(C n - C p ) + j(a n - a p )] x 10 , (17) of 2U8 Pb with the density derivative of the nuclear sym- 



with a = 7.33, /3 = -2.45 fm" 1 , and 7 = -3.62 fin" 1 . 

The parametrization (1171) may be easily understood if 
we consider the PWBA expression of Alr given above 
in Eq. (|16[) . At low-momentum transfer, the form fac- 
tors F n (q) and F p (q) of the neutron and proton densities 
(these arc point densities in PWBA) can be expanded to 
first order in q 2 , so that the numerator inside brackets in 
Eq. $Tjo$ becomes -(q 2 /6) ((r 2 )„ - (r 2 ) p ). In 2pF den- 
sity distributions we have (r 2 ) q — (3/5) C 2 + (7n 2 /5) a 2 . 
Now, assuming constancy of F p (q 2 ) in the nuclear models 
and taking into account that C n + C p ^> C n — C p and 
a n + a p ^> a n — a p , it is reasonable to assume that the 
variation of Alr is dominated by the change of C' n — C p 
and a n — a p as proposed in Eq. (|17l) . 

In the analysis of a measurement of A^r in 208 Pb 
through parametrized Fermi densities, one could set C p 
and a p to those known from experiment [40| and then 



metry energy arises from the bulk part of Ar np . The 
latter result implies that an experimental determination 
of the equivalent sharp radius of the neutron density of 
208 Pb could be as useful for the purpose of constrain- 
ing the density-dependent nuclear symmetry energy as a 
determination of the neutron rms radius. 

We have discussed the shapes of the 2pF distributions 
predicted for 208 Pb by the nuclear mean-field models in 
terms of the so-called "halo-type" (C n — C p = 0) and 
"skin-type" [a n — a p — 0) distributions of frequent use 
in experiment. It turns out that the theoretical models 
can accomodate the halo-type distribution in 208 Pb if 
the symmetry energy is supersoft. However, they do no 
support a purely skin-type distribution in this nucleus, 
even if the model has a largely stiff symmetry energy. 
Let us mention that the information on neutron densities 
from antiprotonic atoms favored the halo-type over the 
skin- type distribution [35, 36|. 
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We have closed our study with a calculation of the 
asymmetry A^r for parity-violating electron scattering 
off 208 Pb in conditions as in the recently run PREX ex- 
periment [2{|, using the equivalent 2pF shapes of the 
models. This has allowed us to find a simple parametriza- 
tion of Alr in terms of the differences C n — C p and 
o-n — Op of the parameters of the nucleon distributions. 
With a measured value of the parity-violating asymme- 
try, it would provide a new correlation between the cen- 
tral radius and the surface diffuseness of the distribution 
of neutrons in 208 Pb, assuming the same properties of 
the proton density known from experiment. 
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